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Abstract 

Computing the topology of an algebraic plane curve means to compute a 
combinatorial graph that is isotopic to ^ and thus represents its topology in R^. 
We prove that, for a polynomial of degree n with coefficients bounded by 2^ , the 
topology of the induced curve can be computed with 0{n^ (n + p^)) bit operations 
deterministically, and with O(n^p^) bit operations with a randomized algorithm 
in expectation. Our analysis improves previous best known complexity bounds 
by a factor of n^. The improvement is based on new techniques to compute and 
refine isolating intervals for the real roots of polynomials, and by the consequent 
amortized analysis of the critical fibers of the algebraic curve. 



1 Introduction 

Problem definition and results We address the problem of topology computation: 
Given an algebraic curve ^(o — V'g_{F) :— {{x,y) e M? \ F{x^y) = 0} implicitly defined as 
the real vanishing set of a bivariate polynomial F e find a planar (straight-line) 

graph G isotopic to -^QThis problem is extensively studied in the context of symbolic 
computation; see related work below. 

We analyze the bit-complexity of the problem. For F of total degree n and integer 
coefficients bounded by 2^ in absolute value, we show that an isotopic graph can be 
computed with 

d{n^{n + pf) 

bit operations with a deterministic algorithm, and with an expected number of 

bit operations with a randomized algorithm, where O means that we ignore logarithmic 
factors in n and p. This is the best known complexity bound for this problem, beating 
the former record by a factor of rp-. 

We give a high-level description of our algorithm first. A more detailed explana- 
tion is given in Section |2] First, V^{F) is transformed to an isotopic Vm(/) by a shear 
such that the sheared curve Vk(./) satisfies certain genericity condition to simplify sub- 
sequent steps. Second, the x-coordinates of critical points are computed as the real 
resultant roots of / and its derivative. Third, the curve is lifted at each critical fiber, 

and G are isotopic if there exists a continuous mapping : [0, 1] x "i^" h-> such that <I>(0, ■) = id^ff, 
Im(<t(l, ■)) = G, and <!>(?, ■) is a homeomorphism between V and Im(<I>(<, ■)) for every t e [0, 1]. 



1 



Figure 1.1: First, the curve is sheared to be located in generic position. Then, critical 
points are projected onto the x-axis defining the x-critical values. In the lifting step, the 
fibers at the critical values and at points in between are computed. Finally, each pair of 
lifted points connected by an arc of C is determined and a corresponding line segment 
is inserted. The right figure shows the final graph that is isotopic to C. 

i.e., the fiber points are computed by real root isolation of the fiber polynomial f\x=a 
(with a a critical x-coordinate), and the index of the (unique) critical point is deter- 
mined. Finally, the number of fiber points in-between critical points is determined by 
the Sturm-Habicht sequence (a.k.a. signed subresultant sequence). The gathered infor- 
mation is sufficient for an isotopic (combinatorial) graph of F; see also Figure [LT] 

Although our algorithm does not differ from related approaches (in fact, its high- 
level description is almost identical to the algorithm described by Gonzalez- Vega and 
Necula 1 16 1), we are still able to derive a complexity bound that improves on all previ- 
ous approaches. This is due to two major novelties in our approach: 

1 . New algorithms for real root isolation of a univariate polynomial ll22l as well as 
for the subsequent refinement of the isolating intervals ||20| express the running 
time in the sum of the local separations and the modulus of the polynomial's 
roots. In particular, when applied to an integer polynomial of degree d and bitsize 
X, the bit complexity for root isolation is bounded by d{d^X^); we use this to 
bound the complexity of the root isolation of the resultant polynomial. Moreover, 
both root isolation and refinement are applicable to arbitrary real polynomials 
by approximating the coefficients and using validated numeric methods. This 
makes them especially useful for computing roots of fiber polynomials, which is 
a critical step in the algorithm. 

2. We consequently use the idea of amortized analysis in this work: When a method 
is applied to each fiber polynomial, we bound the sum of the costs. Usually, that 
sum gives the same complexity as the worst-case bound for a single fiber, which 
means that not all fibers can be bad at the same time. As the main theoretical 
novelty, we bound the complexity of isolating all fiber polynomials by (5(n^p^) 
using this technique. 

Related work Computing the topology of a curve is a problem considered by nu- 
merous papers. We can distinguish existing approaches into such which permit to 
shear the curve as a first step E §1 1.6] lfT2lll23llfT6llfT5]| . and such which do not per- 
mit a shear lfT8]| ll9l. The latter approaches also reveal geometric information about 



the curves, for instance the coordinates of critical points. A mixed approach is taken 
by 1 14 1, where shearing the curve is allowed, but geometric information is still obtained 
by undoing the shear in a post-processing step. 

No matter whether they initially shear or not, almost all mentioned techniques use 
the same three-step approach as the method presented in this paper: they project the 
critical points of the curve, lift the fibers at critical x-coordinates, and connect the fiber 
points by segments corresponding to paths on the curve. The approach by |9| is an 
exception since it avoids computing complete fibers. Instead, it isolates the critical 
points of the curve in and finds an isotopic graph by subsequently subdividing the 
real plane. Another subdivision algorithm by |7| avoids root isolation altogether by 
subdividing regions containing critical points down to a so-called evaluation bound. 
Since root isolation is usually the bottleneck of topology computation in practice, this 
approach looks promising; however, both theoretical and practical comparisons are 
missing to date. 

The time complexity for topology computation has been considered by several of 
the aforementioned approaches. For simpler comparison, let := max{n,p}. Arnon 
and McCallum [ 1 1 gave the first polynomial bound of 0{N^^). Gonazalez-Vega and El 
Kahoui ifTSl improved this to d{N^^) (using classical arithmetic), and Basu et al. ^ 
§1 1.6] prove a bound of 0{N^'^). The best known bound of 0{N^^) has been presented 
first by Diochnos et al. il2l . The same bound was shown also by Kerber ||T9| for the 
algorithm presented in 1(141 . In the two last-mentioned papers, the technique of amor- 
tized analysis is extensively used. In the technical part of our analysis, we sometimes 
refer to those works when the analysis of substeps is identical. 

Outline We start with a more detailed description of our algorithm in Section |2] In 
Section[3] we fix the required notation for the technical magnitudes needed in the com- 
plexity analysis. The analysis starts in Section |4j where we consider univariate poly- 
nomials in general and fiber polynomials in particular, and bound quantities of these 
polynomials like their Mahler measure, coefficient size and separation. In Section |5] 
we summarize the running time of the subalgorithms (e.g., gcd computation and real 
root isolation) needed for the main result. Finally, Section |6] proves the running time 
of our topology algorithm, combining the amortized bounds from Section |4] with the 
subalgorithms from Section [5] 

2 Algorithm 

We start with a description of the topology computation algorithm to be analyzed; 
see Algorithm [T] for pseudo-code. The input is a square-free bivariate polynomial F, 
representing an algebraic curve ^ = yR(F) in M? by its zero-set. The output is an 
embedding of a graph in the plane that is isotopic to 

Initially, the curve F is transformed to / by means of a shear, that is, f{x,y) = 
F{x + sy,y) for some shear factor s e Z. Since the sheared curve C — Vk(/) is isotopic 
to it suffices to compute a graph isotopic to C. The shear factor is chosen such that 
C is in generic position, that is, every critical point has a distinct x-coordinate and there 
exist no infinite arcs that converge to a vertical asymptote. In particular, the leading 
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coefficient of each fiber polynomial f\x=a '■— f{(X,y) G ^[y] at an arbitrary a G C is 
an integer A shear factor satisfying these conditions can be computed deterministi- 
cally or probabiHstically, where the deterministic approach achieves a sHghtly weaker 
asymptotic running time than the probabilistic approach. We refer to the correspond- 
ing paragraph in Section|6]for details. We also remark that this initial step of finding a 
generic shear factor is the only step in the algorithm that is (possibly) randomized. 
In the next step, the x-critical points (i.e., all points p £C with = 0) are projected 

onto the (real) jc-axis via resultant computation. We write /,, :— Let Sres,(/,/J|,) G 
Z [x] the i-th subresultant polynomial and sres, (/, fy) be the i-th subresultant coefficient. 
R :— sreso(/,/y) is the resultant of / and /,. which is the determinant of the Sylvester 
matrix of / and Since / is in generic position, the set VrI^) of real roots ai , . . . , 0^ 
of R contains exactly the projections of all x-critical points. W.l.o.g., we assume that 
«!,...,«,„ are in consecutive order. The set of all points on C located above a certain 
a, is called a critical fiber. The y-coordinates of these points are defined by the roots 
of the critical fiber polynomial /|.r=a, G K^]. 

It is well-known |10| that the curve C is delineable between a, and a,+i, that is, 
C consists of disjoint function graphs which we call arcs from now. Distinct arcs can 
only meet at critical points of the curve, and by genericity, there is exactly one such 
critical point per fiber. By these considerations, the following information is sufficient 
to compute an isotopic graph for /: 

(i) The number of points in each critical fiber, which equals the number of distinct 
real roots of f\x=ai- 

(ii) The index of the unique critical point in each critical fiber, which equals the index 
of the unique multiple root of /|x=a, - 

(iii) The number of arcs between two critical fibers. For arcs in between two critical 
fibers above a,- and o:,+ i, this number equals the number of real roots of fa with 
an arbitrary as («,■, a,+i). 

In all three steps, we use the subresultant coefficients of fiber polynomials. The 
following property |25, §4.4], (|2] §8.3.5] shows that we get them for every fiber by 
evaluating the general subresultant at the corresponding x-coordinate: 

Lemma 1 (Specialization property). For any a € M and any i, 

Sres,(/,/v)|.v=a - Sres,-(/|,-=„,/|;=„). 

For (jij), we first compute the square-free part of each critical fiber polynomial. For 
that, we initially compute the subresultants of / and fy with cofactors, that is, we 
compute M,-,v; S Z[x,y] satisfying 

SreSi(/,/v) = M,/ + v;/j, 

and such that degy(M,) < n — / — 2, deg^,(v,) < n — / — 1. For a critical fiber at a, we 
compute ka '■= dtggcA{f\x=a,f\'x=a) using the well-known property ||2] Prop.4.24] 

ka := mm{k > | sresi(/,/y)(«) ^ 0}. 
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Algorithm 1 Topology computation 
1: procedure Top(F) 

2: Compute 5 S Z such that f{x,y) := F{x~\-sy,y) is in generic position 
3: Compute Sreso(/,/j), . . . ,Sres„(/,/y) with cofactors m;,v; s.t. Sres,(/,/j) = 
Uif + Vify. 

4: Isolate the real roots ai , . . . , of Sreso (/, fy) 
5: for a e {ai,...,am} do 

6: ^^deggcd(/|.,=„,/|;^„) 

7: C v^_i > C|.v=a is the square-free part of f\x=a 

8: Isolate the real roots of C I j=a 

9: Identify the index of the multiple real root of f\x=a 

10: end for 

11: Compute qo,...,qm G Q with qi_\,ai < qt and compute the number of real 

roots of f\x=cij for / = 1 , . . . , in 
12: Construct and return a combinatorial graph isotopic to / 
13: end procedure 



It follows with r?, Prop. 10. 14, Cor. 10. 15] that v<^^_i|v=a is the square free part of 
f\x=a- On this polynomial, we apply the root isolation algorithm from [22 1. The 
results yields the number of real roots and an isolating interval for each root which can 
be further refined to any desired precision. 

For the index of the critical point is computed with the following lemma used 
in IW\- 

Lemma 2. Let Sres, j (/, /,,) denote the coefficient o/Sres, (/,/v)/or y (in particular, 
Sres,-.,(/,/y) = sres, ■(/,/-,,)). For k := ka, define the rational function 

n,. ^ Sl-eSt,A:-l(/,/y)(x) 

Then, the multiple root of f\x=a is Pice). 

Indeed, using this rational expression, /3(a) can be approximated until it can be 
uniquely assigned to one of the isolating intervals of the fiber polynomial. 

Finally, for computing the number of arcs between consecutive critical points ( [mj i, 
we choose rational values qo,...,q,„ with < a, < ^, for all / = l,...,m. The 
number of fiber points at qi can be determined by the signs of sres,(/,/y)(<7,) using 
Sturm-Habicht sequences I.17J . The counting function is easy to compute if the signs 
are known, but its definition is quite lengthy. We refer to lfT4l Sec. 2] for a summary. 

3 Notations 

We fix the following notations and conventions: For a positive real number 0, we write 
Lip : = log ^. We say that an integer polynomial g (uni- or bivariate) is of magnitude 

(c/, A ), if its total degree is bounded by d, and each integer coefficient is bounded by 2^ 
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in its absolute value. For a univariate polynomial g, we denote V{g) the set of distinct 
(complex) roots of R and the multiset of roots of g, that is, each root of g occurs 
as many times in y{g) as its multiplicity as a root of g. 

For a univariate polynomial g € C[x] of magnitude {d,X) with roots zi,.--,Zd, we 
write \cf{g) for the leading coefficient of g. We define the mot bound of g as 

r(g) :=logmax{l,max{|z,| | l,...,d}}, 

the local separation of g at Zi as 

sep(g,z;) := min |zi-Zj|, 

('■;):z,-/z; 

the separation of g as 

sep(g):== min sep(g,z,) 

1=1 a 

(the latter two definition only make sense if g has at least two distinct roots) and 

zev{g) 

The Mahler measure of ^ is defined as 

d 

Mea(g):=|lcf(g)inmax{l,|z,|}. 

1=1 

It is known H Prop. 10.8 and 10.9] that Mea(g) < ||^||2 < Vd+l-2^, and so, 

logMea(^) = 0(A+logt/). (3.1) 
Finally, the local gcd degree of g is defined as 

k{g) :=deg(gcd(^,g'))- 

Throughout the paper, / denotes the bivariate square-free integer polynomial ob- 
tained by shearing our input polynomial F, that is, f{x,y) = F{x + sy,y) with a generic 
shear factor s G Z. The polynomial / is of magnitude (n, t), where n — degF and T 
depends on the bitsize p of the coefficients of F and the bitsize of s. In our analysis, 
we will first compute the bit complexity of our algorithm in terms of the magnitude of 
/ and then relate the result to the magnitude of F. 

As already used in Sectionjij we write := Sres,(/,/,) e 'Z[x,y] for the i-th 
subresultant polynomial and sres, (/,/v) € Z[x] for the i-th subresultant coefficient. For 
convenience, we also write sri := sreSi(/,/v). The resultant polynomial of / and fy 
is defined as R := siq. We can apply Hadamard's bound to immediately read off that 
R is of magnitude (n(n — l),c • n(T + logn)) for some constant c. We further denote 
V{R) := {ai, . . . ,a,-} (with r < n{n — 1)) the set of critical x-coordinates of / and, 
w.l.o.g., we assume that the first m roots ai , . . . , a„, are exactly the real roots of R and 
that they are in consecutive order 

We are mainly interested in the fiber polynomials f\x=a of / with a G C. If a is a 
critical x-coordinate of /, we also talk about critical fibers and critical fiber polynomi- 
als. For shorter notation, we also define 

Fa :=r(/|.^=a), sep„ :=sep(/|.v=a), Ea ■■^T.{f\.x=a), Meaa := Mea(/|.:t=a) and ka :=^(/|x=a). 
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4 Amortized algebraic bounds 

In this section, we investigate the fiber polynomials of / at critical ic-coordinates. For 
various magnitudes, such as root bounds or local separations as defined in|3] we derive 
upper bounds that depend on n and T. We consequently consider all critical fibers at 
once because this leads to the same bounds as considering only the worst fiber among 
the critical fibers. 

Lemma 3 (Mahler bound). Let g be a univariate polynomial of magnitude (t/, A), and 
let V' C ^{g) be any multiset of roots of g. Then, 

^ logmax{l,|a|} <logMea(g) =0(A+logt/). 

In particular, for g — R, the sum is bounded by 0{n{x -\-\ogn)) . 

Proof. Obviously, we can replace V' by ^{g) for an upper bound on the sum. Thus, 

£ logmax{l,|^|}<log n max{l,|^|}<log^^<logMea(g). 



which proves the statement together with ( 3. 1 1. □ 



With this simple result, we can already bound the sum of the root bounds over all 
critical fibers: 

Lemma 4. For any multiset V' C 1^ [R), 

£ logr„ = (9(«2(T + log«)). 

aeV 

Proof. Note that, for any univariate polynomial h — ^4=0 hix', it holds that |25 , Cauchy's 
Bound] 

r(/i) < 1 +max{/!o,...,/!„}, 

so it is enough to bound the coefficients of f\x=a- Notice that every coefficient is given 
by g{<x), where ^ S is a polynomial of magnitude (n, t). It is thus straight-forward 
to see that 

r„ < 1 + («+ l)2"max{l, |a|}" < (« + 2)2^max{l, |a|}", 

and so 

£ \ogTa<n^\og{n + 2)+n^X + n\og ]~[ max{l,|a|}. 

aeV aeV 

The result follows from applying Leinma|3]to the last summand. □ 
Lemma 5. For any multiset V' Q'f (R), 



logMeaa = 0(n^(T + logn)). 



aeV 
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Proof. Notice that Meaa > 1 for every a <E V{R), and that the Mahler measure is mul- 
tipHcative, that means, Mea(g)Mea(/!) = Mea(g/!) for arbitrary univariate polynomials 
g and h. Therefore, 

52logMeaa< £ logMea(/|.,=„) =logMea ]~[ /|.v=a ) ■ 

Considering / as a polynomial in x with coefficients in Z[y], we have that |2, Thm. 4.16] 

-TT ^ res.,{f,R) 

and, thus, 

£ logMea(/|,=„) < logMea(res,-(/,/;)). 

aeV 

It is left to bound degree and bitsize of TeSx{f,R)- Considering the Sylvester matrix of 
/ and R (whose determinant defines res^(/,7?)), we observe that it has n rows with co- 
efficients of R (which are integers of size (9(n(T + logn))) and rows with coefficients 
of / (which are univariate polynomials of magnitude (n, t)). Therefore, the y-degree 
of re Sy (/,/?) is bounded by rr', and its bitsize is bounded by (9(n^(T + logn)). Using 



(3.1 1, this shows that logMea(reSj(/,7?)) = (9(« (r + logn)), as claimed. □ 



Lemma 6 (factorization to multiplicities). R can be decomposed into R^ R\ - •/?„_i 
such that Ri e Z[x] and V (Rt) ^ {a £ V {R) \ ka = i}. 

Proof. Without loss of generality, assume that R is primitive (otherwise, decompose 
its primitive part, and multiply Ri by the content of R). We define Sq :— R, and 5, := 
gcd(5,_i , sr,). By construction, V{Si) = {a E V{R) \ ka > i}. Also, since ka < n for 
all a, degS„ — 0, and thus S„ — I because S„ divides R and R is assumed primitive. We 
define Ri :— It is then straight-forward to verify all claimed properties. □ 

In the subsequent proofs, we require the application of the generalized Davenport- 
Mahler bound that we state here. See ifTSl Thm. 3. 9] for a proof. 

Theorem 7 (generalized Davenport-Mahler bound). Let g E C[t] with n := degg > 2 
and exactly r <n distinct complex roots V := V (g) = {^\ , . . . ,^,-}- J^^t G = {y,E) be a 
directed graph on the roots such that: 

• G is acyclic, 

• for every edge (o;,j3) EE, it holds \a\ < |j3|, and 

• the in-degree of any node is at most 1. 
In this situation. 




min{«,2/i— 2r}/3 



^lc%yfMea(g)'-i 
For the case that G has no edges, the left side simplifies to 1. 
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n 



For the next lemma, recall from Sectionjsjthat L,^ = log ' and sr, = sres,(/, /,.) S 



Lemma 8. For every subset V' C V{R), 



L L,,,^(a) = 0{n\T: + \ogn)) 

aeV 



Proof. We first "complete" the sum by writing 

' I log ' 

aeV(R) 



Next, we show that both summands are bounded by (9(n-'(T + logn)). We first 
prove that LaG\/*log |srA.^(a)| = (9(n^(T + log«)) for any multiset V* C ^(7?). Thus, 
in particular, the second summand in the above equation achieves the latter bound. We 
apply the Davenport-Mahler bound for each f\x=a with a e V* using the empty edge 
set. This yields: 



- V|lcf(/|x=a)|Mea'„"«-' \ma) \^ ) 



min{H,2n— 2ma}/3 



Note that sresj-^ (/!;,.=„, = srA.„(a), that lcf(/|.,=a) = lcfy(/), and that the two 



rightmost factors are both bounded by (i)" from below. Therefore, we have that 



1 > 



V^WT^Mea"-' 
Taking the logarithm of the inverse and summing up, we obtain: 
1 tr- 

1^ E log|sri.„(a)|<-log|lcf,(/)|+« ^ logMea„ +2«3log«. 

The first term is bounded by n^T, and the second term is bounded by 0{n^ (t + logn)) 
by Lemma 5 It remains to prove that Laevffi) log i ,. = Oirv'i.x + logn)). We 
decompose/? = /?i ■ ■ Rn-x according to Lemma |6] and obtain 

■j^ n— 1 

^ '°^R77a^ = ^ ^ -log|srK«)l 

n-1 / \ 



E E -log|sr,-(a)|+ E log|sr,(a)l 

1=1 \aGr(ff,) aer{R,)\V(R,) , 



n-l 

E-log 

1=1 



n srKa) 

aer{Ri) 



n-l 

E E log|srK«)l 

(=1 aenR,)\V{R,) 
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From our above considerations, it follows that the second summand is bounded by 
n^(T+logn)) because the multisets 1 
y{R). For the first summand, we have 



0{rv' [x +\ogn)) because the multisets ^{Ri) arepairwise disjoint and, thus, {Y {Ri)\y {Rt)) C 



n-l 
I 



-log 



n sr^o:) 

aer(Ri) 



n-l 

L 

1=1 

n-l 



log 



res(sr,-,j;,) 

\cf{Ri)'^^&i'-'') 



n-l 



deg(srO log \\cf{R,)\ - log |res(sr,-,7;,)| 



<n^ 
n-l 



>1 



«2log]~[|lcf(/?,)| =n2|lcf(/?)| =0(«^(T + log«)) 



i=l 



□ 



The next lemma bounds the logarithmic inverses of the local separations of an 
arbitrary univariate polynomial. We consider this result to be of independent interest. 

Theorem 9. Let g e M[f] be an arbitrary polynomial of degree d and let k := k{g) = 
deggcd(g,g')- ForV'QV{g), 



Y ^sep(g,,5) 



: 0(t/logMea(g) +L|,sresi(g,g')|)- 



In particular, the bound holds for T.{g) as defined in Section^ 

Proof. Write m := d — k and let , . . . , £,m denote the roots of g. Moreover, let F := 
r(^) > denote the root bound of g. First of all, since every local separation is upper 
bounded by 2^^+', 

2(r+iW]-[sep(^,<^)> n sep(g,^). 

We concentrate on the product on the right hand side first. Observe that, when the i^'s 
are considered as vertices in the complex plane, each sep{g,£,j) is given by the length 
of an edge connecting £,j to its nearest neighbor. This induces a directed graph on the 
vertices, which is known as the nearest neighbor graph ITTl (if a root has more than 
one nearest neighbor, we pick the one with highest index). Let £0 denote the edge set 
of this nearest neighbor graph. We can rewrite: 



n sep(g,<^) 



Our goal is to apply the Davenport-Mahler bound on this product. However, the 
nearest-neighbor graph does not satisfy any of the required properties in general. We 
will transform the edge set £0 into another edge set £3 that satisfies the requirements 
of the Davenport-Mahler theorem, and we will relate the root product of £0 with the 
root product of £3 . 
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Note that a direct property of nearest neighbor graphs is that all cycles have length 2 ifTTI . 
In the first step, we remove one edge of every cycle: 

£i :={(<^„<?,)e£o|/<;V(^„^,-)^£o} 

This removes at most every second edge, and for every removed edge, there is some 
edge in £1 with the same value. Since every root product is bounded by 2""+' from 
above, we can bound 

2(r+i). Yl 1^,-^,1 > n l^.-^'f 

In the next step, we re-direct the edges in £1 in order to satisfy the second condition 
of the Davenport-Mahler bound 

El ■■= {iz„zj) I iizi,zj) e £1 V {zj,z,) e £1) A < |z,| v = |z,| A / < j))} 

In simple words, every edge points to the root with greater absolute value. Note that E2 
does not contain any cycles, because the absolute value of a root is non-decreasing on 
any path, and if it remains the same, the index increases, thus no vertex can be visited 
twice on such a path. Since the only difference between £1 and E2 is the orientation of 
edges, we have 

n n 

Finally, we need to ensure the last condition of the Davenport-Mahler bound, namely 
that each vertex has in-degree at most 1 . For that, if several edges point to some ^j, we 
throw away all of them except the shortest one (in the definition, if the shortest edge is 
not unique, we keep the one with the maximal index): 



Er.^m,^j)eE2\yi^k,^j)eE2:\^k-^j\>\^,^^j\yi\^k-^j\ = \^.-^j\Ak<i)} 

Another basic property of the nearest neighbor graph is that two edges that meet in 
a vertex must form an angle of at least 60° . It follows that the degree of every vertex 
is bounded by 6. Since E2 is a subgraph of the nearest neighbor graph, possibly with 
some edges flipped, the degree of every vertex is still bounded by 6. Since all edges in 
E2 point to the root with greater absolute value, it can be easily seen that the in-degree 
of is even bounded by 3. So, £3 contains at least ^ many edges. Since we always 
keep a smallest edge pointing to a we can bound 

Putting everything together, we have that 

n i^;-^,i>2-5^(r+W n 

{^,4i)£Eo \ii„ii)eE, 
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Zi3 meets all prerequisites of the Davenport-Mahler bound and we can thus bound 



nsep(g,^)=2-^(r+^) Yi i^,_^,i>2-«^(r+') n l^.-^'i 



#£3 ^ ^ 



^lc%^Mea(g)'"-i \m \mj \^ 



> 2 



-6rf(r+i) 



V|sres,(g,gO| /I 
y|lc%yfMea(g)'' ■ \d 



(in the last term, we have simplified some of the factors which are irrelevant for our 
argument). Passing to the inverse and taking logarithms, we obtain 

L ^sep(g,,5) < 6(i(r+l)+3L|5resi(^y)|+31oglcf(g)+6t/logMea(g) + 12t/logJ 

= 0{d\ogMea{g) +L|sresttey)| +t/r + loglcf(^) +t/log<i), 

and the last three terms are all dominated by dlogMea{g), because Mea(g) is larger 
than 2^, lcf(^), and logii by definition. □ 

Let y' C V{R) be the set of all roots of R where f\x=a has at least two roots. It has 
been shown lfT3l Prop. 3.73] that 

L Lsep„ =0(«^(T + l0g«)). 

aeV 

We will prove that this is also true when replacing sep„ by the (strictly larger) E^. 

Theorem 10. Let V' C V{R) be the set of all roots of R, where f\x=a has at least two 
roots. Then, 

E„ = (9(«^(T + logn)). 

aeV 

Proof. For fixed a € V' , we denote nia :^n — ka > 2 the number of distinct roots of 
f\x=a- We apply Theorem|9]on each f\x=a (all of degree n) to obtain 

Y^^a = 0(«logMea„+L|srj,j„)|) = C>(n ^ logMeaa+ ^ L\s.v,^(a)\) 
aeV aeV aeV aeV 

The first sum is bounded by (9(n^(T + logn)) (Lemma|5]l and the second sum by 
0(n^(T + logn)) (Lemmalsll. □ 



5 Basic algorithms 

In the whole paper, all complexity bounds for algorithms refer to the bit complexity, that 
is, the number of bit operations needed to achieve the algorithmic task. Our bounds 
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usually depend on the magnitude (n,T) of the input polynomial. For simplicity, we 
mostly ignore logarithmic factors in n and T in the complexity bounds, and write O to 
refer to bounds where logarithmic factors are omitted. We assume asymptotically fast 
multiplication on integers, hence, multiplication of two «-bit integers has a complexity 
of ()(«). 

Basic operations We list the complexity of several basic operations on uni- and bi- 
variate polynomials next. We omit most of the proofs - see fl9', §2], f2, §8] and ll24l 
§11.2] for a more complete treatment. The maybe most fundamental non-trivial sub- 
operation that we need in our algorithm is the evaluation at rational values. 

Lemma 11 (Rational evaluation). f/ ITPl Lemma 2.4.10]) Given g e l^x] of magnitude 
{d,X), and a rational value ^ such that c and d have a bitsize of at most O. Then, 
evaluating g{^) has a complexity of 

d{d{X+da)). 

Another fundamental operation is to compute the greatest common divisor of uni- 
variate polynomials. 

Lemma 12 (gcd computation). Let g,h e Z[x] be both of magnitude {d,X). Computing 
their gcd has a complexity of 

d{d^X), 

and the resulting gcd has degree at most d and its coefficients have a bitsize of 0{d + 
A). 

Closely related to computing the gcd is the square-free part of a univariate polyno- 
mial, which is given by g/ gcd{g,g'). 

Lemma 13 (square-free part). Let g G Z[x] be of magnitude {d,X). Its square-free part 
g* can be computed in 

d{d^X) 

and it has degree at most d. The bitsize of each coefficient of g* is bounded by 0{d + X). 

Root isolation Given a univariate polynomial, we want to compute its real roots. 
By "computing", we understand to compute a list of isolating intervals, each interval 
containing exactly one root of polynomial. For this subtask, we use the result from ll22l 

Theorem 14 (Root isolation). Let g = Y!i=i8i^' ^ U^H square-free polynomial 
with \g„\ > 1, r := T{g) the root bound of g and E :— E(^). Then, we can compute 
isolating intervals for the real roots of g in time 

(5(t/(c/r+E)2). 

For that, every coefficient must be approximated to a precision of 

(5(t/r+E) 

bits after the binary point. 
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However, isolating intervals are not always sufficient for our algorithm; we often 
need that, in addition, each interval is smaller than a given e > 0. In this context, 1201 
study the problem of root refinement: 

Theorem 15 (Root refinement). With the same notation as in Theorem^M^ assume that 
isolating intervals for the real roots of g are known, and let £ > be an arbitrary real 
value. Then, computing isolating intervals of g of width at most £ needs at most 

did{dr + I.)^+d^U) 

bit operations, and each coefficient must be approximated up to a precision of 

OiLe+dF + I.) 

bits after the binary point. 

Putting both results together, we obtain 



Theorem 16 (Strong root isolation). With the same notations as in Theorem 14 given 
a polynomial g and £ > Q, we can compute isolating intervals of g of width at most £ 
within at most 

d{d{dT + I.f+d'^U) 
bit operations, and each coefficient must be approximated up to a precision of 

d{u+dT+i.) 

bits after the binary point. 

The special case of integer polynomial has been considered in the aforementioned 
papers, too. A bound of 

d{d^X^+d^Le) (5.1) 

has been shown for this problem. Being slightly more careful, we obtain the same 
bound also for non-square-free polynomials: 

Theorem 17 ((Strong) root isolation, integer case). Given a polynomial g €E not 
necessarily square-free, of magnitude {d,X), we can compute isolating intervals for 
the roots of g with at most 

d{d^X^) 

bit operations. If the intervals are additionally required to be of width at most £, they 
can be computed with a number of bit operations bounded by 

d{d^X^+d^Le). 

Proof Let g* denote the square-free part of g. By Lemma [T3] it can be computed 
within 0{d^X) bit operations and its magnitude is {d,d + X). Using (5.1 1 forg* would 



yield a worse complexity than claimed. Instead, we use the bounds from Theorem 14 
and Theorem 16 Note that k{g*) = deggcd{g* ,{g*)') = and so. Theorem |9] yields 
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e 0(t/logMea(^*)+L|si.es|j(g* = 0{d\ogMea{g*)), where the last equality 

follows from sreso(g*, (g*)') > 1 because g* and its derivative are integer polynomials. 
Moreover, Mea{g*) < Mea{g) because g* divides g over the integers. It follows that 

e 0{d{X +log«)) — d{dX). Moreover, because g and ^* have the same roots, we 
can apply the Cauchy bound on g to get r(^*) — T{g) = d{X). Plugging in everything 



in Theorem 14 and Theorem 16 yields the desired bounds. □ 



In some situations, we do not require small isolating intervals, but rather the con- 
trary: we seek for rational values which separate the roots of the polynomial from each 
other and have a small accumulated bitsize. The following result achieves this; its proof 
is a direct consequence of the properties of the isolating intervals returned by the root 
isolation algorithm from ll22l . 

Theorem 18 (Intermediate values). For an integer polynomial g of magnitude {d,X) 
with m real roots zi , . . . we can compute rational values qo, . . . ,qm with < 
Zi < qi and bitsizes yQi-.-^ym that sum up to 0{d{X + \ogd)), performing not more 
than d{d^X^) bit operations. 



Proof. The algorithm from 11221 uses classical bisection to compute isolating intervals 
{ai,bi) for the real roots zi of g with 

^^^^<|«,_^,.|<2Jsep(^,zO, 

see II22I Theorem 18]. Thus, the bitsize of the endpoints of a, and bj is bounded 
by logL(y(g j..) +log(16t/^). For qi := ^'~2^"' , the bitsize y of qt is also bounded by 
0(Lo-(g.z,) +logii). Thus, summing up Yi over all / yields an upper bound of (9(E(g) + 
dlogd) = 0{d{X+logd)). □ 

Note that, in particular, ,^,) is an isolating interval for z,-. 

Interval arithmetic The main operation that we will perform on an algebraic number 
is: Given h E Z[x], a G M algebraic and 5 > 0, compute some r E Q such that \r — 
h{a)\ < 5. In other words, we want to approximate h{a) to absolute precision Lg. 

We achieve this task by using interval arithmetic: For two intervals Ii — [ai,bi], 
h = [aiM]^ we set 



*B(/l-/2) 



= [a\+a2,b\+b2] 
= [01-^2,^1-02] 

= [mm{aia2,bia2,aib2,bib2},vaax{aia2,bia2,aib2,bib2}] 



^ih/h) :- «(/i-[^,^]),if0^/2. 

02 bl 

For a polynomial h = LJLq fl,x'^ and an interval /, we evaluate according to the Horner 
schem^ 

<B(/z(/)):=«8(flo+/-(fli +/•(...))) 



^It should be noted that, unlike in 1201 . we use exact interval arithmetic, that is, the boundaries are not 
rounded to a floating point grid. 
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where each a, is interpreted as interval [a,-, a,]. We observe that *B (/!(/)) contains the 
image of h under /, although it can be much larger than that. Also, note that an el- 
ementary arithmetic operation in interval arithmetic consists of at most 4 elementary 
operations on the interval boundaries; therefore, we can still use asymptotically fast 
methods for interval arithmetic. In particular, if the boundaries of the interval are ra- 
tionals with bitsizes bounded by a, we can evaluate 58(/z(/)) with d{d{X +da)) as in 



Lemma 1 1 (with h being of magnitude {d,X)). 



Going back to the problem of approximating h{a) to precision Lg, assume that a 
is given by some isolating interval / of size e (initially set to j). We evaluate h{I) 
using interval arithmetic to obtain an interval J — ^h{I) which contains h{a). If the 
diameter of J is smaller than 5, any value in the interval yields a valid approximation 
value. Otherwise, £ is set to and the method is repeated. 

To quantify when / becomes "small enough", we use a technical result on interval 
arithmetic. 

Lemma 19. ( I]19\ Lemma 2.5.20]) Let h g "Llx] be of magnitude (t/,A) and I be an 
interval of width < £ < 2. Then, for each a € I and each y € *B/z(/), we have 

\y-h{a)\ <2''£2^max{l,|a|}^"'. 

Theorem 20. Let g,h Cz be of magnitude (t/, A). Let ai , . . . , o:„, be the real roots 
ofg, 5i , . . . , 5,„ e M such that < 5; < 1, and 5 :— Y[T=i ^i- Then, approximating h^ai) 
to precision Sjfor all i = 1, . . . ,m has a total complexity of 

d{d^X^ + d^Ls). 

Proof. Let /; be the isolating interval of a, . If /; is refined to size 



£/ : = 



5i 



2d+^2^ max{l,|a,|}''-i 
the distance of y € *B /(/,) to /i(a,) is bounded by 

\y-Kad\<\5, 



using Lemma 19 and by the triangle inequality, the length of ^h{I) is smaller than 5,. 

Thus, /,■ must be refined at most to precision £,. Note that 5,- > 5 for all /, thus it 
suffices to refine each /, to size 

e:=2- 



'2''+i2^max{l,|a,|}''-i' 
Since |a,| e 0{X), we can bound 

Li;^0{Ls+d + X+dX)=d{Ls+dX) 
Refining all /, 's to size £ takes 

did^X^+d^L^)) = did^X^+d^Ls)) 
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bit operations by Theorem 17 which is the desired bound. 

It is left to argue why the interval evaluation and the failing tries with too large 
values of e in the algorithm do not increase the complexity. Note first that if strong 
root isolation is applied for the same polynomial and decreasing values of e, the cost is 
determined by the call with smallest e in the sequence. Furthermore, since e is squared 
in every step, the bitsizes of the interval boundaries are doubled in each iteration. Thus, 
the evaluations are essentially determined by the last evaluation, where the boundaries 
have abitsize of 0{X +L£.). Therefore, the final evaluation step for a,- costs 0{d{dL£. + 
A)). We show that 



Y,U,^0{Ls+d{l+\ogd)). 

1=0 



Indeed, 



£ Le, = £ L5, + (^Z + 1 ) log 2 + A log 2 + (^/ - 1 ) log n max{ 1 , 1 a,- 1 } , 

(=0 (=0 (=1 

and the latter is bounded by 0{X +logii) by Lemmajs] Thus, the interval evaluations 
are bounded by d{d^Lg +d^X), which is dominated by the overall complexity bound. 

□ 

The previous proof can be used for a slightly more general result: We will not just 
approximate /i(a, ) for a single h, but for a whole sequence /i, 1, . . . all of same 
magnitude. Instead of just multiplying the above bound by k, we can do better: 

Theorem 21. Let g, ai , . . . , a,„, 5i , . . . , 5„, and 5 be defined as before. Moreover, let 
e^'.yOtl m denote a set of m ■ k polynomials, all of magnitude {d,X). Then, approxi- 
mating hi j{(Xi) to precision 5ifor all i — I, . . . ,m and j — I, . . . ,k has a total complexity 
of 

d{d^X^+k{d^X+d^Ls)). 

Proof. The previous proof shows that, once a, is refined to precision e,-, the width 
of *B/j(a) is less than or equal to 5, for any h of magnitude {d,X). Thus, we still 
need not more than 0{d^X^ + d'^Lg) bit operations for the refinements, no matter 
how many /i/./s we consider. The additional summand k{d^X +d^Ls) arises because 
we have to bound the cost of the interval evaluations: We have shown the bound of 
0{d{dLe. +X))) for evaluating a single /i/.y; since there are k polynomials to evaluate, 
the evaluations cost are 0{kd{dL£. + X))) for a,-. The results follows from bounding 
the sum of Lg. in analogy to Theorem [20] □ 

6 Topology Computation 

Theorem 22 (Main result). Algorithm^has a bit complexity of 

d{n\n + pf) 

if the shear factor is chosen deterministically, and an expected bit complexity of 
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if the shear factor is chosen probabilistically. 



Generic position We first ensure that the sheared curve V{f ) =V{F{x^- sy,y)) is in 
generic position, that means 

• deg(/) = deg^.(/) (the leading coefficient of /, considered as a polynomial in y, 
is a real value) 

• for each a G M, f\x=a has at most one multiple root. 

Geometrically, this is equivalent to the absence of vertical asymptotes and covertical 
critical points. Original and sheared curve are known to be isotopic, so that computing 
the topology of the sheared curve is sufficient. 

We follow the approach from ifTSl Lemma 16]: For a bivariate polynomial /, we 
denote A(/) as the discriminant of the square-free part of the resultant of / with /, . If 
a curve V{f) is generic, A(/) 7^ must hold. Thus, the goal is to find a i G Z such that 

A{F{x + sy,y))^Q. 

Theorem 23. f/72) Lemma 16]) Given a polynomial F of total degree n and integer 
coefficients ofbitsize p, deterministically computing an s £7^ such that f{x,y) = F(x + 
sy,y) is in generic position has a bit complexity of 0{n^{n + p)). 

The proof computes A for an indeterminate s using a fast subresultant algorithm ETl . 
and evaluates A(/(x + sy,y)) G 1,[s] at 0{n'^) integer values using multipoint evalua- 
tion lIM] Cor. 10.8] to identify a non-root of A. 

Alternatively, we can also compute a generic direction with a probabilistic method 
(which is also recommended in practice for efficiency reasons). Since the number 
of "bad" shear factors (such that yield a non-generic curve) is bounded by n'* + n lfT9l 
p. 1 17-1 18], we can simply take an integer i from the range {1, . . . ,2(«^+«)} at random 
and check whether A{f{x + sy,y)) is zero or not for this specific s. If it is zero, another 
s is chosen. Note that at least half of the values in the range must be non-roots, so the 
expected number of iterations is 2. 

Computing A(F) requires (5(«^(n + p)) bit operations, because computing the re- 
sultant requires (5(n^p) operations [21] and the resultant is of magnitude (n^,n(p 



log«)). Computing its square-free part needs (5(n^p) bit operations with Lemma 13 
and the square-free part is of magnitude [n^ ^n{p -\- n)) . Computing its resultant, and 
thus its discriminant requires d{n'^ {n + p)) bit operations 1.19. Thm.2.4.16]. We can 
conclude 

Tlieorem 24. Probabilistically computing an s such that f(x + sy,y) is in generic 
position has an expected bit complexity of OirP {n + p)). 

The resulting shear factor s is ofbitsize O(logn), thus f{x + sy,y) has a maximal 
coefficient size of T = 0{p +log«) = 0{p). From now, we assume that / has been 
transformed into generic positions in all subsequent steps. In particular, the results 
from Section]?] apply for /. 
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Computing subresultants and critical values The computation of the subresuhant 
polynomials Sreso(/,/y), ... ,Sres„ (/,/,) with their cofactors can be done in 0{n'^T) 
bit operations ['2 1 j . Each Sres,- (/, /v ) is a polynomial of x-degree at most n^, y-degree at 
most n — /, and maximal coefficient size of n(T + logn). In particular, R := Sreso(/,/y) 
is a univariate polynomial of degree n^, and its roots are the critical x-coordinates of 
the curve. Computing them is now an application of Theorem 17 to which yields a 
complexity of 

d{n^ ■ (n(T + log«))2) ^ d{nh^) = d{n^p^). 



Computing the A;'s Recall that for a root a of R, we denote by ka the degree of 

gcd(/|,=„,/|U). 

Theorem 25. Computing all ka 's has a bit complexity of 

Proof, ka is defined by the minimal index k such that sres4.(/,/v)(a) 7^ 0. Checking 
whether STesk{f,fy){a) vanishes can be done by computing gcd(/?, sres^^ (/,/,)), and 
checking whether the sign of the gcd changes when evaluated at the boundaries of any 
isolating interval for a. Since both polynomial are of degree (at most), and their 
coefficient bitsizes are bounded by n(n + T) (n(logn + T) for sreSi.(/,/-y)), one such 

We need to do this 
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gcd operation has a bit complexity of 0{n [n + t)) by Lemma 
at most n times. 

18 to choose the evaluation points. Let m<n^ denote the number 



We use Theorem 



of real roots of R, and let q'o , . . . , q,n denote the rational intermediate values. Computing 
them requires (5(«^T^) bit operations. Let 70, ... , Jm denote the corresponding bitsizes. 
We have to evaluate each gcd at each value qj. One such evaluation costs 

d{t?{n{n + T)+t?Yj)), 
and the total costs ai^e therefore bounded by 

m m 

6{nY^n^(n{n + x)+n^yj)) = d{n^{n + x)+n^ ^ Jj). 
j=o ' j=o 

Because the 7,'s sum up to d{n^t), we obtain a bound of d{n^t) for this step. □ 



Computing the fibers We next bound the costs for isolating the roots of the fiber 
polynomials. 

Theorem 26. Given f, R, and CCi, . . . ^CCr as above. Assuming that the square-free part 
gi of f\x=aj is known for / = 1, . . . ,r, isolating the real roots of all of them is bounded 
by 

d{n^X^)^d{n^p^). 

Proof. We have to show two parts. On the one hand, we have to bound the running time 
of the root isolation algorithm, assuming that a sufficient precision of the coefficients 
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is available. On the other hand, we need to bound the time for computing a sufficient 
precision. 

In the proof, we will write f*\x=a for the square-free part of f\x=a- Note that 
/* \x=a = Ci\x=a (for some /) where C; G "^[x^y] is a cofactor polynomial of a subresul- 
tant of / and /,,; see Prop. 10. 14, Cor.10.15]. Let C, j G denote the coefficient 
of Ci at yK It is known that each C, j is a polynomial in x with degree at most n^, and 
bitsize at most n(T + logn). 

For the first part, recall from Theorem[T4]that the running time of root isolation for 
f*\x=a is 

d{n{nT{f*\x=a)+nf%=af)- 

We observe that r(/*|,=a) = Ta and T.{f*\x=a) = J^a- Moreover, with Theorem|9] we 
have Ed G (9(nlogMeaa + .sr<;^ (a)). Thus, we obtain a bit complexity of 




(5(«^(nr„, +nlogMea„, +L,,^^ 
1=1 ' 

= o[n'(^Ya}j +«^|^tMea„,^ +ni^L, 

The first sum is dominated by the second because Ya^ < Meaa.. The second sum 
is bounded by C>(n^(T + logn)) according Lemma [s] The last sum is bounded by 
C>(n^(T + log«)) by Lemma^ After all, we get a complexity of (5(n^T^) for this step. 

For the second part, we use the second part of Theorem [l4] Let 5, such that Lg. 
is the number of bits to which the coefficients of f\x=aj need to be approximated for 
isolation. Because L5. = 0{nYaj +Ec(;), it can be seen by the same methods as above 
that the L5 's sum up to 6>(n^(T + logn)). Moreover, let C, be the cofactor polynomial 
of / and fy that defines the square-free part. Our problem is to find approximations of 



Q.o(o:,), ■ ■ .jCi^nif^i) with a precision of 5,. We can use Theorem 21 to bound the costs, 
setting hi j ^ Cij, d i~ n^, X ^ n(T + logn), k ^ n and L5 -s— n (r + logn), which 
yields 

for getting sufficient precision for root isolation. □ 



Detecting the multiple root Let a:—ai be a critical x-coordinate, and jSa.i , . . . ,l5a.mi 
the roots of f\x=a- Since / is in generic position, exactly one of the jSaj 's is a multiple 
root. Since we have worked with the square-free part of f\x=a in the isolation, we do 
not know yet which root is multiple. Recall from Lemma [2] that the multiple root is 
given by j3 (a) with 

s ^ sres^.^_i(/,/v)(x) 
k-siesk,k{f,fy){x)' 

We describe a simple algorithm to find the index of the multiple root: We set £ := j and 
refine / until 7 := *Bj3 (/) has a width of at most £. We also refine the isolating intervals 
of the fiber polynomial f\x=a to size e. If J overlaps with only one isolating interval of 



20 



f\x=a, we have found the multiple root. If there is more than one such overlap, we set 
e to and retry. 

It is not difficult to see that the above algorithm terminates at the latest when e < 
^sep„. We next prove a bound on the width of / such that this is guaranteed. For 
simpler notation, we set pa sresA.^,A.^_i (/,/,.) e Z[x] and qa ^"^^^ka.kaif Jy) G 
1[x]. 

Lemma 27. If the width of I is smaller than 

\q{a)\he^a 



25+r„ (2«'2«(^+i°g«)max{l, |a|}«')^ ' 



the width of ~ — -^^jj^ is smaller than ^sep„. 

Proof. Note that p and q are of magnitude (n^,n(T + logn)). Let / be isolating for a 
with width smaller than 5a- Set y e ^p{I). By Lemma [T9{ we have that 

y-pia)<e:= '^^ ^' ^" 



25+ra 2«('=^+'°g") max{ 1 , 1 a I }« 

and the analogous inequality holds for ^(a). 
£ has the following three properties: 

1 . e < is^^ . Indeed, we can rewrite e as 

\q{a)\ 1 sep„ \q{a)\ 



E 



2 8 2ra+i 2"'2"(-^+i°g")max{l,|a|}«" 
and the latter factors are both smaller than 1 . 



2. e < 



\ti{a^sep„ ^ same argument as injlj and noting that Fa > 0. 



3. e < ^'^^2\p{a)\" • ^S^i'^' replace Ta by and exploit that 

\pia)\ < 2"'2''(^+'°s"Vax{l,|a|}"' 



Fix some y E *8j3 (/). We can write y as 



y = 



Pia)+ei 
k{q{a)+e2) 



with |ei|, \e2\ < £■ So we get that 

1 p{a) pia)+ei 



\P{a)~y\ 



< 



q{a) q{a)+e2 
e2p{a) 



< 



e2p{a) 



q{a){q{a)+e2) q{a)+e2 



q{a){q{a)+e2] 



q{a)+e2 



< 



£\pia)\ 



I?(«)ll(?(«)+e2)| \qia)+e2\ 
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By[l] we have that |^(a) +e2| > ^^M^ thus 



\q{a)^\ \q{a)\ - 16 16 8 

using [2] and [3] Thus, it follows by triangle inequality that two values in Sj3(/) cannot 
have a distance of more than . □ 

Lemma 28. Let V' C V{R) denote the real roots ofR. Then, 

£ Ls,^0{n\t + n)). 

aeV 

Proof. Note that < 5a < 1 for any a G C. Thus, we can bound 

aeV aeV{R) 

< 5+ E r„+2n4+2n3(T + logn)+2n2 ^ max{l,|a|} + 

aeV{R) aeV(R) 
1 



aeV(R) a£V(R)^^'<a\('^) 



The first sum is bounded by 0{n^{z + logn)) (Lemma|4|, the second sum by (9(n(T + 
logn)) (Lemma I3J, the third by 0(«^(t + logn)) (Theorem 10 1, and the fourth by 
C»(«^(T + logn)) (Lemma[8]l. □ 

Theorem 29. Identifying the multiple roots for all fibers can be done in 

dinh^) = din^p^). 

Proof. Recall the algorithm to find the critical point. It consists of three major building 
blocks: refining the isolating interval / of a (to size da in the worst case), evaluat- 
ing ^8)3(7) using interval arithmetic, and refining the isolating intervals of the fiber 
polynomials to a size of ^sep^, in the worst case. 

We analyze each part separately. For the first part, it is enough to refine each 



isolating interval of 7? to a precision of Y-aev' ■ Using Theorem 17 and Lemma 28 
this can be done with at most 

d{nh^+n'^{T + n))^d{nh^-) 

bit operations. For the second part (interval arithmetic), we note that the costs are dom- 
inated by the last evaluation because e is squared in every iteration. The bitsizes of the 
interval boundaries are bounded by 5a, so that the last evaluation has a bit complexity 
of 

6{n^{n{T + logn)+n^LsJ. 
Summing up over all a's yields 0{n''{n + t)). Finally, we bound the third part (re- 



fining the fiber polynomials) using Theorem 16 Note that, with the notation of that 

Theorem, L^ i ^j^, and isep^ is dominated by La which also appears in the bound. 

It follows that the term "d^Lg" is dominated by the first summand, and the complexity 
reduces to the cost of isolating the fiber polynomial, which is bounded by (5(n^T^) with 
Theorem |26l □ 
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Fiber points at intermediate positions The last missing step is to compute the 
number of arcs between two ascending critical x-coordinates. We do so by com- 
puting the number of fiber points over a rational x-value between these critical x- 



coordinates. Recall from Theorem 18 that we can find such rational values qo,. . . ,q„j 
in time (5(«^T^) = d{rfip^), and their bitsizes sum up to d{rr'x). 

The number of roots of the polynomial f\x=qj is determined by the sign pattern of 
the principal subresultants of f\x=qj and its derivative according to the Sturm-Habicht 
sequence |17|. Let ji denote the bitsize of qj. Evaluating the n principal subresultants 
at qi has a cost of 

d{tr'{n{'Z + \ogn)+n-yi)), 



by Lemma 11 and summing up over all qi yields a bit complexity of 0{n x) = 0{n p). 

To summarize, we have shown that every step in Algorithmjljis bounded by d{rfip^) 
(except for deterministically computing a shear factor), which finally proves our main 
Theorem I22I 



7 Conclusion 

Our work has proven a new worst-case bound for the reference problem of computing 
the topology of an algebraic curve. The result would not have been possible with- 
out improving the complexity of real root isolation f22l and root approximation EOl : 
however, we emphasize that none of the algorithm that achieved the previously best 
complexity bounds |[T2llfT9l would improve to our bound if just the real root isola- 
tion and refinement algorithms are exchanged. This shows that the careful amortized 
analysis performed in our work is an integral ingredient for the obtained result. 

A natural question would be how to further improve the result. To get substantially 
lower bounds than the presented one, we believe that deeper insights on the algebraic 
properties of algebraic curves are necessary. For instance, a bottleneck in the current 
analysis is the isolation of the resultant polynomial which is assumed to be a general 
polynomial of magnitude («^,«t). However, a counting argument on the dimensions 
shows that not every polynomial of that magnitude can appear as the resultant of a 
curve of magnitude (n,T), which leads to the question: is the isolation of a resultant 
polynomial possibly easier than for a general polynomial? At the same time, it might 
be worth to think about lower bounds on the problem of topology computation; to our 
knowledge, no lower bound except the trivial n(n^) (complexity of a planar graph with 
vertices) is known. 

One might also ask about the practical quality of the presented algorithm. Note 
that our algorithm is very similar to the AlciX algorithm 1 19|| 14| which has been im- 
plemented as part of the algebraic kernel package of CGAlj^|4|; the main difference 
is the root isolation at fiber polynomials: while our methods computed the square-free 
part of the polynomial for isolation, AlciX avoids this computation by using the m-k- 
Bitstream Descartes which is a variant of the Descartes method that can cope with one 
multiple root in the fiber The reason for this choice was better practical performance 



The Computational Geometry Algorithms Library, http: //www. cgal . org 
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compared to the computation and isolation of the square-free part, so we do not expect 
our method to be faster than AlciX in practice. Moreover, 1 6 1 and 1 3 1 have recently pre- 
sented new approaches which generally outperform AlciX. It is an interesting question 
whether the same complexity result as in this work can be achieved for AlciX, or for 
one of the two most recent methods. 

A related, but less studied question is the complexity analysis for computing the 
triangulation of an algebraic surface. An algorithm for this problem has been presented 
by 1 5 1, where computing the topology of the projected silhouette curve is a crucial 
building block. Since that curve is of magnitude («^,«t) (for a surface of magnitude 
(n,T)), a complexity bound of (5(n'**T^) appears possible, and we pose the question 
whether this bound can really be achieved. 
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